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DEFINITION OF SYMBOLS 


Symbol 

Definition 

A 

Vector potential 

a 

Radius of cylinder 

B 

Magnetic induction 

c 

Velocity of light 

D 

Electric displacement 

E 

/■v/ 

Electric field 

H 

Magnetic field 

J 

Current density 

k 

Propagation constant (2-rrn/X) 

L 

Energy flux density at a surface 

C 

Length of cylinder 

A A AAA 

£, m, n, r, e 

Unit vectors 

M 

Magnetization 

m 

Refractive index 

P 

<•*«/ 

Polarization vector 

P 

Dipole moment 

S 

Poynting vector 

s 

Scattering matrix 

Si 

Scattering amplitude functions 

V 

Velocity 


v 



DEFINITION OF SYMBOLS (Concluded) 


Symbol 

a 

a 

0 

e 

eo 

d 

X 

M 

Mo 

P 

a 

da 

d£2 

<P 

<t> 

CP 


Definition 

Polarizability 
Tilt angle 
Bisectrix angle 

Permittivity or dielectric constant 

Permittivity of free space 

Scattering angle 

Wavelength 

Permeability 

Permeability of free space 
Charge density 
Electric conductivity 

Differential scattering cross section 
Scalar electric potential 

Azimuth angle between the plane of incident beam and z-axis and 
the scattering plane 

Angular frequency 
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ELECTROMAGNETIC SCATTERING BY CYLINDERS - AN INTRODUCTION 

SUMMARY 


This report, the first of a series of papers, presents an introduction to the 
scattering of electromagnetic waves by cylinders. The basic scattering features particular 
to cylinders are discussed from the Rayleigh-Gans theory. 

In the first part of the report a review of electromagnetism and a derivation of 
dipole radiation are presented. In the second part, the basic scattering theory is given. 
From the equations for dipole radiation, the theory of scattering by cylinders according 
to the Rayleigh-Gans theory is developed. 

A computer program to calculate the scattering intensity functions is presented. 


INTRODUCTION 


This report is the first of a series of papers describing the scattering of electro- 
magnetic waves by cylinders. The applications are numerous, e.g., in studying cylindrical 
plasma, polymers, paints, rodlets, and platelets; in the fields of astronomy, chemistry, 
meteorology, and physics. Cylindrical scattering differs from spherical scattering in 
several important ways. There is a loss in the degree of symmetry introducing anisotropic 
processes which can generally be characterized by a shape factor. There are off-diagonal 
components in the scattering matrix which arise from the mixing of the electric field 
components. 

The general methods to be discussed in these reports are (1) Rayleigh, (2) 
Rayleigh-Gans, (3) infinite cylinders, (4) finite cy Under s by neglecting end effects, (5) 
geometric scattering, (6) Watson transformations, (7) Wiener-Hopf technique, (8) per- 
turbation techniques, (9) symmetry techniques, (10) S-matrix theory, and (11) diffraction 
theory. 



ELECTROMAGNETIC THEORY 


Maxwell Equations 

Light, electromagnetic radiation, is described by Maxwell’s equations. All 
theoretical optics problems, including scattering problems, can be solved formally by 
using Maxwell’s equations. 

The state of excitation (force) which is established in space and detected by the 
presence of an electric charge is said to constitute an electromagnetic field. The electro- 
magnetic field is a vector field, and associated with it are the two vectors E and B, 

/v rsu 

called the electric field vector and the magnetic induction vector, respectively. The 
spatial and temporal derivatives of these two vectors are defined by Maxwell’s simultaneous 
partial differential equations [ 1 ] : 



( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


where p is the electric charge density. The vectors D, H, and J , are the electric 

displacement, the magnetic vector, and the electric current density, respectively. For 
macroscopic media, the dynamical response of the aggregates of the atoms is summarized 
in the constitutive relations for isotropic, permeable, conducting dielectrics [2] : 


D = 

eE 

J = 

ctE 


B 

H = 



Ai 


( 5 ) 

( 6 ) 
( 7 ) 


2 





where e is the electric permittivity, a is the specific electrical conductivity, and p is 
the magnetic permeability. In this work, e and p are assumed to be scalars. 


Maxwell’s equations are the results of particular experimental laws: equation (1) 
from Coulomb’s law, equation (4) from Ampere’s law, equation (3) from Faraday’s 
law, and equation (2) from the observed nonexistence of magnetic monopoles. 

Table 1 gives the units of the electromagnetic quantities. 


TABLE 1. ELECTROMAGNETIC QUANTITIES IN SI UNITS 


Symbol 

Quantity 

Units 2 

E 

Electric Field 

N/C 

H 

Magnetic Field 

C/m-sec 

D 

Electric Displacement 

C/m 2 

B 

Magnetic Induction 

N-sec/C-m 

S 

Poynting’s Vector 

N-m/sec-m 2 

c 

Speed of Light 

m/sec 

e 

Electric Permittivity 

C 2 /N-m 2 

M 

Magnetic Permeability 

N-sec 2 /C 2 

a. N = newtons, 

C = coulombs, and m = meters. 



Wave Equations 

Maxwell’s equations predict the existence of electromagnetic waves propagating 
with the velocity of light, which lead to the electromagnetic theory of light. The wave 
equations are obtained by vector operations [3] on Maxwell’s equations. The equations 
are usually solved in applications using the appropriate boundary conditions for a 
particular problem. Assuming the isotropic case (the constitutive equations), 

B 

H = — 

~ M 

D = e E , and J = a E 


3 











V X (VX E) = V(V- E)- V 2 E 
/■«*/ 

but V-E = 0 in the charge free region. Therefore, equation (13) becomes 



Similarly, take the curl of equation (11): 

VX (VX B) = ff(iVX E +iie (VX E) (15) 

~ ~ 9t ~ 
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Then, upon substituting from equation (10): 

9B 9 2 B 

V X (VX B) = -an -ne — f 

~ 9 1 9 r 

or 


(16) 

Equations (14) and (16) are the inhomogeneous vector wave equations for E and B. 

In a nonconducting medium, a = 0 and the second term in equations (14) and 
(16) vanishes, leaving a propagation equation for waves with a velocity v = l/V/Te . 

The boundary conditions at an interface between two media are [4] : 

1 . The normal component of B is continuous. 

/v 

(B 2 - Bj ) • n = 0 

2. There is a discontinuity in the normal component D equivalent to K, the 
surface charge density, 

(D 2 - Dj ) • n = K 

3. The tangential component of E is continuous, 

(E 2 - E, ) X n = 0 

^ iV> 

4. There is a discontinuity in the tangential component of H equal to L, the 
surface current density, 

(H 2 -H,) X n = L 
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Plane Waves 


We now consider the homogeneous vector wave equations (a = 0) in which the 
field depends only on one coordinate r and on time. Such waves are said to be “plane”: 

a 2 E 

V 2 E -pe = 0 

~ at 2 

a 2 B 

V 2 B - jue ^ = 0 

~ at 2 

Assume 

E (r , t) = E 0 f(£, t) 

where f(r , t) satisfies the scalar wave equation: 

V 2 f (r , t) - p e f (r , t) = 0 . (17) 

Now let 

f (r , t) = f (n • r - vt) ; 

then 

Vf (r , t) = f' (n • r - vt) n • Vr 
but Vr is the idemfactor [3] . Thus, 

A __ A 

n -V r = n 
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and 


V f (r , t) = n f (n • r - vt) 

V 2 f (r , t) = V • V f = V • n f = n • V f = n • n f" = f" 

Also 

f=v 2 f" (18) 

From (17) and (18) we have 

v = -J- (19) 

Hence, the solution for E and B is 
E (r , t) = E 0 f(n-r - vt) 

/“V/ 

B(r,t) = B 0 g (n • £ - vt) 

The wave equations do not relate f, g, E, and B, but Maxwell’s equations do. 

From Equation (8) 

V • E = 0 


7 



and substituting 


8 


( 20 ) 


;tion of 


propagation. Similarly, for B 


( 21 ) 





3> l W 


we have 


-B — - B 0 g - - B 0 — - — 
~ Vue 

hence, 

-B 0 g' 

— = n X E 0 f' 

Vjue" 

— B 0 dg = nXE 0 df 
s/ye ~ 


B - Vmc n X E 


( 22 ) 


and B are mutually orthogonal and are perpendicular to the direction of propagation 


Propagation Constant 

For a plane wave propagating along the positive z-axis, the component of E 
along the x-axis can be written as: 


E x = Aei( kz - Wt ) 


(23) 


where only the real part of the exponential is to be taken as the physical quantity. 
Substituting equation (23) into equation (14) yields 


-k 2 


A e^ w ^ " kz )+ga ^ico A e^ w ^ ' kz ^ + ge c u 2 


Ae i(cot-kz) = o 


9 




Therefore, the propagation constant, k, is related to the constitutive constants of the 
medium 


k 2 = peco 2 +ijuaco 


and may be represented as [4] 


k = a x +i0i 


where 


a x = co 



( 




and 


J3i = co 


1/2 




( 24 ) 


Thus, for a conducting or dissipative medium (a > 0), the propagation constant is 
complex, providing a damping part e'^ 1 z to the wave which corresponds to absorption. 


Polarization [2] 


The plane wave 
E = e x E 0 f (n r -vt) 


is a wave with its electric field vector always in the direction e x . Such a wave is said to 

be linearly polarized in the direction of ei . To describe a general state of polarization, 
two linearly polarized independent waves are needed. Consider harmonic plane waves 

i(k • r - co t) 
f (r, t) = e ~ 


10 




where 


= v , and k = kf 
k ~ 

Then the general state of polarization can be described with the following two waves: 


ik • r - ioj t 
E x = e, E t e~ ~ 

ik • r - ico t 

E 2 = e 2 E 2 e 


with 


Bj = yfixe n X Ej j = 1,2 


and perpendicular to e 2 . The amplitudes E, and E 2 are complex numbers to 

allow the possibility of a phase difference between the waves. The general harmonic 

wave solution is a linear combination of Ej and E 2 , 

/*>»/ 

i(k • r - cj t) 

E (r , t) = & Ej + e 2 E 2 ) e ~ ~ (25) 


If E x and E 2 have the same phase, equation (25) represents a linearly polarized wave, 
with its polarization vector making an angle 6 - tan -1 (E 2 /E t ) with c, and having a 
magnitude E = \/e, 2 + E 2 2 . 

If Ej and E 2 differ by a phase of 90 deg, then 


i(k • r - cot) 

E (r,t) = E 0 (^ ± i e 2 ) e ~ ~ 


and the wave is circularly polarized. The “+” sign is for left circularly polarized waves 
or is said to have positive helicity. The negative sign is for right circularly polarized 
waves or negative helicity. 
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For polychromatic radiation, the harmonic plane wave solution has the form 


E (r , t) 



e n (k) 


i(k • r - w t) 

IV /v 

e 


( 26 ) 


where E 0 „ are amplitude functions determined by E (r , t) at a specified time, and 
e n is the polarization vector. 


Poynting Vector 

For a periodic field such as the electromagnetic wave, the energy crossing a unit 
area per unit time is given by the Poynting vector 


S = E X H 


The time-average flux of energy is given by the real part of the complex Poynting 
vector: 


S = 1 (E X H*) 


Because 


( 27 ) 


B 

H = — = nX E 

~ At A* ~ 

by equation (22) for plane waves, then 



The energy density £ is 
IS I 

£ = 

c 
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The Poynting vector is in the direction of propagation and can be regarded as a 
measure of the intensity of the wave at a point. The Poynting vector gives the flux 
density at the surface, L, which is called by Van de Hulst the intensity, “I”, 


S=c£k=kL=k “I” 


Electromagnetic Potentials [2] 

Vector potentials are often introduced to simplify problem solving. According 
to equation (2), 


V - B = 0 


then the field of B is always solenoidal. Since 


V ■ (VX v) = 0 


B can be represented as the curl of another vector A 0 : 


B = V X A 0 


(28) 


However, A 0 is not uniquely defined by equation (28) because B is equal also to the 
curl of some vector A , 


B = VX A 


(29) 


where 


A = A 0 - V 
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Since 


V X (Vtfr) = 0 

i p is any scalar function of position. If B is replaced in 
9B 

V X E = - — 

at 

by either equations (18) or (29), then 

vx|e + 

or 

vx (g- 

Thus the fields of the vector 

3A 0 
E + 

~ at 

or 

9A 

E + 

at 

are irrotational and equal to gradients of scalar functions <t> and 0 O : 
3A 0 

E = - V0 O - 


9t / ° 
3A \ 

— = 0 
at / 


(31) 
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and 


E = 


- V 0 - 


9A 

~af 


(32) 


The functions 0 and 0 O are related by 


0-0 o + 


90 

9t 


To show this, recall that A is given by 


A = A 0 - V 0 


and from equations (31) and (32), 


9A 0 9A 9A 0 

V0 O + — = V0+ = V0+ - V 

at at at at 


V 0o — V 0 - V 


90 

at 


or 


V 0 


V0 O +V 


90 

at 


and finally, 


0-0 o + 


90 

at 
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The functions A are vector potentials of the field, and the 0’s are scalar 

'V 

potentials. Using D = e E and B = p H, then 


D = - el V 0 + 


3A 


at 


H = ~ VXA 
~ M 


Equations (1) and (4) of Maxwell’s equations yield: 

a 2 A 


V X V X A + jue +ue — 

- at at 2 


M J 


V 2 0 + V 


3A 


at 


Imposing the Lorentz condition, 


V ■ A + ixe = 0 

- at 


we then have 


(33) 


V 2 0 - jue 


a 2 ^ 

at 2 


= v • Ao + IX e 


30o 

3e 


(34) 


where 0 O and A 0 are particular solutions. 0 and A are defined by 


3 2 A 

VXVXA-W-A + ix6 — =ixl 

- - at 2 


V 2 0 - — f = 

at 2 


(35) 
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Now, using V X V X A = VV • A - V • V A , 



( 36 ) 


For the homogeneous case 


D = - V X A' 


E = - — V X A' 
~ e ~ 


H = -V*' 


3 A' 

~aT 


B = -At V 0' + 


3 A 


at 


the wave equations for these potentials are: 


3 2 A' 


V 2 A' - pe 


at 2 


- 0 


(37) 


v a 0 '-„ e = 0 


(38) 


with the condition that 


V-A’+Ate = o 

at 


Retarded Potentials 

Starting from equations (35) and (36) 

3 2 A 


V 2 A - Ate 


at 2 


= - At J 


(39) 


(40) 
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V 2 0- M e = - -i- 
3t 2 e 


P 


(41) 


a solution of the inhomogeneous linear equations for an initial value problem can be 
represented as the sum of the solution of these equations without the right-hand side 
(RHS) and a particular solution of these equations with the RHS. To find the particular 
solution, we divide the whole space into infinitely small regions and determine the field 
produced by the changes located in one of these volume elements. Because of the 
linearity of the field equations, the actual field will be the sum of the fields produced by 
all such elements, i.e., an integral [5] . 

The charge de in a given volume element is, generally, a function of time. If we 
choose the origin of coordinates as the center of the volume element, then the charge 
density is p = de(t) 5 (R) where R is the distance from the origin and 8 represents 

the Dirac delta functions. Thus, we must solve the equation 

v2 0-pe-^l = d e(t)8(R) (42) 

3r e ~ 

Everywhere 5 (R) = 0 we have 
V 2 0 -pe = o 

at 2 u 

For the case of central symmetry, i.e., 0 is a function of R only, we have 



Assuming 0 = x(R, t)/R; then 

/v 


#X_ 

3R 2 



= 0 
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But this is the equation for plane waves, whose solution has the form 


x = t ‘(*-f) +f ’( t+ f) ■ 

where v = 1 Isfile . Since we want only a particular solution of the equation, it is 
sufficient to choose only one of the functions f t and f 2 . Usually it is convenient to 
take f 2 = 0, since this physically represents a wave that is observed before it is generated. 
Then, everywhere except at the origin, <p has the form, 


0 = 



(43) 


So far, the function x is arbitrary; we now choose it so that we also obtain the 
correct value for the potential at the origin. We must select x s° that at the origin 
equation (42) is satisfied. This is easily done, noting that as R 0 the potential increases 
to infinity, and, therefore, its derivatives with respect to the coordinates increase more 
rapidly than its time derivative. Consequently as R-> Owe can, in equation (42), neglect 
pe3 2 0/3 1 2 compared with V 2 <p. Then equation (42) goes over to the familiar equation 
leading to Coulomb’s law, i.e., from 

V 2 0 = — — de(t) 5(R) 
e ~ 

and since 
we have 


0 = -Ml) 

4?reR 
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Thus, near the origin equation (43) must go over into the Coulomb’s law, from which it 
follows that x(t) = de(t)/4rr e , that is, 



For an arbitrary distribution of charges p(x, y, z, t), one can find the solution 
to equation (41). Let de = p dV and integrate over the whole space. To this solution 
of the inhomogeneous equation we can still add the solution <p 0 of the homogeneous 
equation. Thus, the general solution has the form: 

<Hx, y, z, t) = -J— - J — p(x', y', z', t - ^ dV'+0 o (44) 

v 

where R 2 = (x - x') 2 + (y - y') 2 + (z - z') 2 and dV' = dx' dy' dz'. R is the distance 
from the volume element dV' to the “field point” at which we determine the potential. 
Similarly, 

f j(x',y',z',t-^-W 

A(x, y, z) = -f- / ^ W + a 0 (45) 

~ 4-n J R ~ 

V 

Without 0o and A 0 the above are the “retarded potentials.” The volume integrals 
in equations (44) and (45) represent the contributions from charge distributions contained 
in the volume element. A 0 and 0 O represent the contributions from fields whose origin 
is external to the system. 

The solution of the wave equation can be accomplished “directly” with Green’s 
functions, G, where G satisfies [5] 


V 2 G - pe = 5(R)5(t) 

and the solution to this equation can be obtained by Fourier transforms. 
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Dipole Radiation [2] 


In the following, the radiation from an electric dipole is discussed. Assume that 
the distribution of charges changes little during the time r'/c where r' refers to the position 
of the source point. Let the radiation of the system have periods of order T. Let a be the 
order of magnitude of the system. The time r'/c is of the order a/c. So that the distribu- 
tion of charges in the system shall not change significantly during this time, it is necessary 
that a/c«T; but cT = X, the wavelength, hence, a condition for the following discussion 
is 


a « X 


(46) 


Assume that 


P (r , t) = p (r) e'* 60 1 (47) 

J (r , t) = J (r) e" iwt (48) 


and that the electromagnetic potentials and fields have the same time dependence. Since 
the time-dependence factor is the same in all terms, it can be dropped. Then from 
equations (45) and (48) the retarded vector potential is 


A (r) 




ik Lr - r' I 
e ~ ~ 


d 3 r' 


(49) 


where lr - r'l = R, k = co/c, and d 3 r' = infinitesimal volume element. Confining the 

source to a region that is small compared with the wavelength and considering r » a 
(Fig. 1), then 


lr-r'l 


r - n • r 


+ 
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SOURCE POINT 



Figure 1. Geometry for the source point, field point, and origin. 

n is a unit vector in the direction of r , and 

-ik (n • r' + . . .) 
e ~ d 3 r’ 

n • r' 

1 - + .. . 

r 



( 50 ) 


If r » a and a « X we can expand the exponential: 


,-ikn • r' + . 


1-HH- + ... 


1 + 


- ik^ (n • r') + 




The m™ term of the vector potential becomes 


A m = 


iae ikr 

4nr 


(-ik) 

m! 


m 


1 + 


ai 

ikr 


a m 

(ikr) 


-)/ 


J(r')(n -r') m d 3 r' 


where aj are integers. In the far field where r » a we have 


lim A m -> 
kr->^ 


pe 


ikr 


47rr 


(-ik) m 

m! 



J(r') (n • r’) m 


d 3 r' 


( 51 ) 
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Taking m = 0, equation (51) becomes 

A(r) = f J(r')d>r' 

Integrating by parts, the integral becomes 

f d 3 r' J (r') = - fr' (V' • J ) d 3 r' = -iw f r' p (r') d 3 r' 

I /■>*» rv J in# rw J rv/ ~ 

The last step comes from the continuity equation 

. JA. = v • J 

at 

or 

iccp = V- J 

/V 

Hence, the vector potential is 


A(r) ' - i "“ eikr P 
~ ~ 4ttt ~ 

(52) 

where p the electric dipole moment is defined as 


P = f £ p (T) d 3 r' 

(53) 


The electric dipole fields are 


B - V X A = ~^ k - (n X p) / 1 - — j— ) 

~ ~ 4n ~ r \ lkr J 
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E = — VXB = -i ^VXB 

~ wjxe ~ k ~ 


_ k 2 
4ne 


(nXp)Xn 

<N# 


e ikr 

r 


+ 


1 

4ne . 


3 n (n • p) - p 


ik \ ikr 

IT) 


In the far-field zone, 



(54) 


Scattering Amplitude 

The scattering of plane waves by a particle is completely described by a knowledge 
of the fields at every point in space. Generally, the observer is interested in knowing the 
fields at his location. He desires to know the scattered fields at large distances from the 
scatterer where the scattered waves are less complex and appear spherical. In this region 
the amplitude of the scattered waves decreases inversely with the distance r from the 

particle, and the phase has a simple e + * kr dependence. The scattering amplitude defines 
the strength and phase of this outgoing spherical wave. Since the amplitude and phase of 
the scattered radiation differ in various directions from the scatterer, the scattering 
amplitude js / is a function of position. Since the electric field is transverse at large dis- 
tances, we write the asymptotic form in the two-component notation: 


ik • r i^i- 

E (r) — * E 0 e (0,*) (55) 

~ ~ r~>oo ~ j- 
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The scattering amplitude is most conveniently described by a scattering matrix 
which allows the polarization directions to be handled in a simple manner. 

Figure 2 shows the coordinates used to describe the direction of the incident and 
scattered radiation along with its polarizations. The x, y, and z axes are fixed to the 
scatterer but chosen so that the z axis is along the direction of the incident radiation. 
The unit vector n is in the direction of the incident radiation having the propagation 


/ 

/ 




k 1 




vector k, and the unit vector n' is similarly defined for the scattered radiation having 
the propagation vector k' [6]. The vectors n and n' (or k and k') define the 

'V /N/ 

scattering plane, and the angle between n and n' is the scattering angle B (0 < 9 < 1 80). 

The components of the electric field are referred to the r - 2 system for the 
incident radiation and the r' - 2' system for the scattered radiation. These systems 
simplify the discussion of polarization direction. The unit vector r is perpendicular to 
the scattering plane, and the x-axis of the particle makes an angle 0 with 2. The unit 
vector 2, which is in the scattering plane, is chosen so that r X £ = n, and the unit 

vector 2' is chosen so that r'X 2' = n\ The unit vectors r and r' are the same. The 
quantities E 0 f and E 0 g are the amplitudes of the components of the electric field in 

the incident beam along the two directions defined by ? and 2. Likewise, E r and Eg 
are the electric field components of the scattered radiation along r' and 2'. 


The scattering amplitude tensor which relates the incident and scattered waves is 
defined by 


j^sca '-S^(B , 0) 


gikr 

r 


- S(0,0)E o 


e +ikr 

r 


(56) 


or 






e ikr 

r 


The four scattering amplitude functions S t , S 2 , S 3 , and S 4 are all functions of 6 and 
0 and, in general, are complex. As stated above, the components (Eg, E r ) and 

(E 0 g,E 0r ) are referred to different sets of coordinates, i.e., (2', r') and (2, r). Hence 
the scattering amplitude tensor is a function of k and k', i.e., 


S(fl,0) = S (k, k') 

— — rss 


The linear relation imphed by the S-matrix (scattering matrix) between E and E 0 
is a consequence of the linearity of Maxwell’s equation. Because E and E 0 are vectors, 
S (k, k') must be a tensor as shown. For no scattering, S reduces to a unit matrix. 
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Differential Cross Section 


The differential scattering cross section da/dQ. is the ratio of the flux scattered 
from the object per unit solid angle to the incident flux on the object per unit area. 


da _ energy scattered/unit time/unit solid angle 
dO energy incident/unit time/unit area 


In terms of Poynting’s vector, S, for the incident and scattered fluxes, we have 


da 

dO 


q 

£sca 

q. 

pnc 


R 2 


k sca 

k inc 


R 2 


Scattering by Electric Dipoles (Rayleigh Scattering) 

We now consider an electric dipole at the origin of a Cartesian coordinate system 
upon which a monochromatic plane wave impinges from along the negative z-axis 
(Fig. 2). 

If the electric dipole is generated by electrons and if the incident wave is a 
harmonic plane wave of the form, 


E = E 0 ie ikz ‘ iwt 


then the equation of motion is 


mr = e E (z = 0) 


Since p = e r , the electric dipole moment is found from 

/v /v 

P = — E (z = 0) , 

~ m ~ 

to be 


oj 2 m 


E 0 i 
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Note that the use of E at z = 0 in the equation of motion implies that r « A, i.e., 
the dipole approximation. 

In a dielectric in the presence of an externally applied electric field, the molecular 
dipoles are streteched and become oriented. A net electric moment is produced in the 
direction of the applied field. The dipole moment per unit volume for an isotropic 
dielectric medium is related to the external electric field by 


P 


= a E 


where a [C 2 m/N] is the proportionality constant, the polarizability, and is dependent 
on the frequency of E. From the discussion on dipole radiation the radiation fields for 

an electric dipole are [eq. (54)] 


E = - 


co 

k 


B X n = 


e ikr 

4neT 


k 2 ^ (p X n) X n ) 


B 


= co pke ikr 

47rr 


w A 

p X n 


where E, B, and p were assumed to have a harmonic time dependence. 
Since the incident electric field is along the z-axis 


line 


= E 0 e 


ikz - ico t 


A 

i 


or 

Sine * (E„ r J + E 0 ,i)e ikz -i“* 
where E 0r = E 0 sin <j> and E 0 g = E 0 cos0. Hence 
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B c „ 0 = - C0M - kelkr (a E inc X n') 


"sea 


4nr 


where we have used 


£ a Sine 

I z - 0 

Then 

B sca = - ae ikr (E 0r fX n' + Eo C ^X n') 

Now since 

A A; A.. 

r = r r X n = - £ 

£ X n' = r cos 0 £ X n = r 


where we have used \ = L (£, n') and cos 0 = sin £. Therefore, 
Ssca = - (E 0 r (-£') + E 0 g (r cos 0 ) 

and 


(58) 


"sea 


to 

k 


Ssca * n 


ae ikr k 2 
47rer 


E 0r (- n' X £') + E 0 g cos 0 (n' X r') 


) 


. at iltr f 
47rer 


(Eo r r' + E 0 g cos 0 £') 


(59) 
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Hence, 



where 


s 2 = 

a k 2 cos 6 


4ne 

s, = 

a k 2 

4ne 


( 60 ) 


( 61 ) 


The time-averaged Poynting vector for the scattered wave is 
^sca “ ^Ssca x JJsca) 

In the far radiation zone, E and H are perpendicular to each other and to the direction 

/v 

of propagation. For a harmonic wave under these conditions, 

H* = H- n X E* ; 

~ V M 

thus 

^sca ~ ~2 (^sca x ( n x J^sca)) J — , 
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or 


_ 1 


;sca 


a 2 k 4 


M (47re) 2 r 2 


(E 0r n' + E 0 gcos 2 9 n') 


(62) 


Therefore, the energy flux densities L = 1/2 \fefy. E 2 for the incoming and scattered waves 
are related by 


a 2 k 4 


"sea 


(4-ne) 2 r 


2 ,2 v inc 


(Line + cos2 e L inCe ) 


(63) 


If <j> = 0 deg, then 


= — - — ^ - cos 2 6 
sca (47re) 2 r 2 mc 


if 0 = 90 deg, 


t = a 2 k 4 L . 
L sca (4ne y r 2 me 


and if 0 = 45 deg or the incident radiation is unpolarized, 


L 


sca 


a 2 k 4 
(47re) 2 r 2 


-J (1 +COS 2 9 ) L inc 


(64) 


For natural incident radiation, this is the same as averaging over all polarization angles, 


L 


sca 


a 2 k 4 
(47re) 2 r 2 



sin 2 0d0 + Lj nc 


cos 2 6 
2ir 



a 2 k 4 
(4 7r e ) 2 r 2 


Y (1 +cos 2 0)L inc 
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This scattering is called Rayleigh scattering. The differential scattering cross section becomes 


(65) 


Rayleigh-Gans Scattering [7] 

Assume that each volume element in some finite volume gives Rayleigh scattering 
and does so independently of the other volume elements. The waves scattered in a given 
direction by all these elements interfere because of the different positions of the volume 
elements. To calculate the interference effects we have to refer the phases of all scattered 
waves to a common origin of coordinates and then add the complex amplitudes (Fig. 3). 




Figure 3. Scattering by a finite volume. 

Consider the volume element in Figure 3. If p is the dipole moment per unit 
volume, the magnetic induction for the scattered field at point P is 



(pdVXm) . 
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Integrating over the volume, this expression becomes 


B = - J*JL f dV ( p x m) (66) 

~ 4tt J R ~ 

With n, a unit vector in the incident direction, and m, a unit vector in the scattered 
direction, the phase of e^^ (Fig. 4), i.e., kR is k ^R 0 + r • (ni - n)J . The wave “2” 

has to go (r • n) farther to get to the scatterer, but for wave “1” the scatterer is (r • ni) 
farther in front of the origin, R 0 is large compared with the body size. Hence, 


B = (p X m) f- \ 
~ ~ \ 4 it / 


e ikR 0 f lk r • (m - n) 

^ — / e dV 

Ro J 


which differs from Rayleigh scattering by a factor of the form JV 5 where 
5 = k r • (m - n). We now evaluate the normalized factor 



Figure 4. Phase change due to scattering. 
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r («,*) = y y eis dv 


(67) 


for a cylindrically shaped body. As seen from Figure 5 the vector m - n has the length 
2 sin (0 / 2 ) along the bisectrix of the direction m and - n. Hence, 


5 = s 2 k sin — 
2 


where 


s = r • 


(m - n) 


Im-nl 


( 68 ) 


(69) 


Hence, 


R (0 ,0) = 


V 



s sm 


0 

ds 


(70) 


where we are integrating over “slices” perpendicular to the bisectrix, and each slice has an 
area D and thickness ds. 
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We have assumed that the phase shifts of the wavelets passing through the medium 
are small, i.e., 


2 ka Im - 1 1 « 1 


(71) 


where a is a length of the order of the size of the particle, and m = the index of 
refraction. 


Circular Cylinders of Finite Length 

The function R(0 , 0) can be calculated for the circular cylinder of finite 
length. Let its length by £ and its diameter 2a, and let the phase shift be small for a ray 
traversing the cylinder in any direction. The orientation of the cylinder with respect to the 
incident wave is arbitrary (Fig. 6). 


INCIDENT 

BEAM 



Figure 6. Geometry of scattering by a cylinder. 
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The volume integration will be performed by using circular slices which are perpen- 
dicular to the cylinder axis since this is easier than using slices perpendicular to the 
bisectrix. 

ikr-(m-n) 

To calculate e let r = i x + j y. Fix the axis of the cylinder so that the 

bisectrix lies in the y - z plane. Then, m-n = a I j + a 2 k and 


r • (m - n) = y a t 

Since Im - nl = 2 sin 0/2 and is in the same direction as the bisectrix. 



or 

a 

aj = 2 sin — - sin p 
2 

where (3 is the angle between the cylinder axis and the bisectrix, and 


r • (m - n) = 2 y sin sin (3 


With this relation and with the phase referred to the center of the disk, 


r d (0 


1 f ik ~ ' 

= D J ' 


(m-n) 


dD 
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becomes 


Rj) (0,0) 


D 



a 

2iky sin — sin (3 
e 2 dx dy 


2 _ 

D 



a 

2iky sin — sin f? 
e 2 Va 2 - y 2 dy 


( 72 ) 


We can make a change of variables by letting 


y = aw 


and 


v = a2k sin — sin (3 
2 


then 


Rp (0 ,0) 


2 

D 



e ivco a 2 yj \ - u* dto 


14- I e ivw doj 

^ 3 J-l 




r 1 


_2_ 

7T 

J 

! e* vw -s/l - dco - 

/ e lvco Vl -«' 2 d«' 



0 

•A) 
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Let 03 ' = -gj in the second integral 



and hence 

Rj) (0 >0) 



= F(v) = -1- J,(v) 
v 


where J n (v) is a Bessel function of order n and 



( 73 ) 


( 74 ) 


The phases of the disk have been referred to the center of the disk. In considering 
the whole cylinder, the phase can be referred to the center of the cylinder. This introduces 
another phase factor given by 

kz k • (m - n) 
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But since 


A A A A 

(m - n) = aj j + a 2 k 


the phase is k za 2 where a 2 = cos ( 3 . This phase effect is given by 
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The final phase factor is the product of the two phase factors, i.e.. 


R(M) 


F ( 


2 ka sin 


sin /3 


M 


kC sin 



The angle p and be related to the angle a between the incident light and the z-axis 
the angle 0/2 and the azimuth angle <j> between the incident light — z-axis plane and the 
incident light-scattered light plane by (Fig. 7) 


cos P = - cos a sin 


e_ 

2 


+ sin a cos 


0 _ 

2 


cos 4 S 


(77) 


AXIS OF 
THE CYLINDER 


INCIDENT 

DIRECTION 



Figure 7. Scattering angles. 




To find the scattering amplitude tensor we follow the steps as in the Rayleigh scattering 
section to obtain: 


B = - (p x j - ) e 1 k R ° v r (0,0) 

47T ~ Ro 


E = upL ( (p x m) X m ) VR (0,0) 

~ 4t / R 0 


( 78 ) 


= k - (aV - } e ^ Ro - (E 0 r r' + E 0 g cos 0 £') R (0 , 0) 

4tf 6(j Rn 


Hence, 



(79) 


We note that since R (0 , 0) is independent of polarization effects the polarization will be 
that of a Rayleigh scattering particle: 


P(0) = 


ISi I 2 - IS 2 I 2 
IS! I 2 - IS 2 P 


1 - cos 2 6 
1 + cos 2 0 


sin 2 6 
1 + cos 2 0 


The differential scattering cross section is 


(80) 


da 

d£2 


(aV ) 2 k 4 
(4tt e 0 ) 2 


R 2 (0,0) -L (1 + cos 2 0) 


(81) 


The appendix gives a computer program listing for calculating the functions F and 
E and the scattering intensity functions. Also included, but not derived in this paper, is 
a phase factor for spheres G. Table 2 tabulates some values of these factors. Figure 8 
shows them for comparison. 
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TABLE 2. THE DIFFRACTION FUNCTIONS 


u 

e<u> 

.00 

1.00000 

.10 

.99833 

.20 

.99335 

.30 

.98507 

.40 

.97355 

.50 

.95885 

.60 

.94107 

.70 

.92031 

.80 

.89670 

.90 

.87036 

1.00 

.84147 

1.10 

.81019 

1.20 

.77670 

1.30 

.74120 

1.40 

.70389 

1.50 

.66500 

1.60 

.62473 

1.70 

.58333 

1.80 

.54103 

1.90 

.49805 

2.00 

.45465 

2.10 

.41105 

2.20 

.36750 

2.30 

.32422 

2.40 

.28144 

2.50 

.23939 

2.60 

.19827 

2.70 

.15829 

2.80 

.11964 

2.90 

.08250 

3.00 

.04704 

3.10 

.01341 

3.20 

-.01824 

3.30 

-.04780 

3.40 

-.07516 

3.50 

10022 

3.60 

-. 12292 

3.70 

-. 14320 

3.80 

-. 16102 

3.90 

-. 17635 

4.00 

-. 18920 

4 . 1 C 

-. 19958 

4.20 

-.20752 

4.30 

-.21306 

4.40 

-.21627 

4.50 

-.21723 

4.60 

-.21602 

4.70 

-.21275 

4.80 

-.20753 

4.90 

-.20050 


F(U» G<U) 

.00000 1.00000 

.99875 .99900 

.99501 .99601 

.98879 .99103 

.98013 .98409 

.96907 .97522 

.95567 .96446 

.93999 .95185 

.92211 .93745 

.90211 .92131 

.88010 .90351 

.85619 .88411 

.83048 .86321 

.80311 .84089 

.77421 .81723 

.74392 .79235 

.71237 .76633 

.67972 .73928 

.64613 .71132 

.61174 .68256 

.57672 .65310 

.54123 .62306 

.50542 .59256 

.46945 .56172 

.43349 .53064 

.39768 .49946 

.36217 .46827 

.32711 .43719 

.29265 .40632 

.25892 .37579 

.22604 .34568 

.19414 .31609 

.16334 .28712 

.13374 .25886 

.10543 .23139 

.07850 .20479 

.05304 .17913 

.02910 .15447 

.00675 .13088 

-. 01397 .10840 

-.03302 .08708 

-.05038 .06697 

-.06602 .04809 

-.07995 .03046 

-.09217 .0141 1 

-.10269 -.00095 

-.11154 -.01473 

-.11876 -.02721 

-.12437 -.03842 

-.12845 -.04836 



U 

5.00 

5.10 

5.20 

5.30 

5.40 

5.50 

5.60 

5.70 

5.80 

5.90 

6.00 

6.10 

6.20 

6.30 

6.40 

6.50 

6.60 

6.70 

6.80 

6.90 

7.00 

7.10 

7.20 

7.30 

7.40 

7.50 
7.60 

7.70 

7.80 

7.90 

8.00 

8.10 

8.20 

8.30 

8.40 

8.50 
8*60 

8.70 

8.80 

8.90 
9.00 
9.10 
9.20 

9.30 

9.40 

9.50 
9.60 

9.70 
9.80 

9.90 


E(U) 

-.19178 
-.18153 
-. 16990 
-. 157 D 3 
-.14310 
-. 12828 
-.11273 
-.09661 
-.08010 
-.06337 
-.04657 
-.02986 
-.01340 
.00267 
.01821 
.03310 
.04720 
.06043 
.07266 
.08383 
. 09386 
. 10267 
.11023 
.11650 
.12145 
. 12507 
. 12736 
.12833 
. 12802 
.12645 
. 12367 
.11974 
.11472 
. 10870 
. 10174 
.09394 
.08540 
.07620 
.06647 
.05629 
.04579 
.03507 
.02423 
.01338 
. 00264 
-.00791 
-.01816 
-.02802 
-.03740 
-.04622 


F(U) 

-.13103 
-.13219 
-.13201 
-.13055 
-.12791 
-. 12416 
-.11940 
-.11374 
-.10725 
-.10005 
-.09223 
-.08389 
-.07513 
-.06606 
-.05676 
-.04734 
-.03787 
-.02846 
-.01918 
-.01012 
-.00134 
.00709 
. 01509 
.02 26 2 
.02963 
.03607 
.04190 
.04709 
.05163 
.05549 
.05866 
.06114 
.06293 
.06403 
.06447 
.06426 
.06343 
.06200 
.06002 
.05751 
.05451 
.05108 
.04726 
.04310 
. 03865 
.03395 
.02907 
.02405 
.01895 
.01381 


G(U) 

-.05705 
-.06453 
-.07083 
-.07598 
-.08002 
-.08300 
-.08498 
-.08599 
-.0861 1 
-.08539 
-.08389 
-.08168 
-.07882 
-.07537 
-.07141 
-.06699 
-.06219 
-.05 70 7 
-.05169 
-.0461 2 
-.04 04 1 
-.03463 
-.02883 
-.02 30 6 
-.01737 
-.01 182 
-.00644 
-.00127 
.00 36 5 
.00829 
.01262 
.01661 
.02025 
.02 35-2 
.02640 
.02890 
.03099 
.03269 
.03400 
.03491 
.03544 
.03560 
.03541 
.03488 
.03403 
.03288 
.03146 
.02979 
.02790 
.02580 
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CYLINDRICAL SCATTERING FUNCTIONS 





Figure 8. Cylindrical functions. 



Randomly Oriented Cylinders [7] 


To obtain the differential scattering cross section in the Rayleigh-Gans domain for 
randomly oriented cylinders, an integration over the solid angle of the orientation of the 
cylinder is performed : 


R 2 


4n 


/ R 2 (0,tf>) 


dfi 


(82) 


where R 2 is an average cylindrical function for random orientation. The values 0 and <p 
are considered fixed. The solid-angle integration reduces to an integration over (3. For the 
case of thin rods (R ->■ E), 


E 2 (z cos 0) d cos p 



_1 

z 



(83) 


where z = k 8 sin 1/2 9. 


The first term is the sine integral. For the case of disks, 


R 2 = 



F(z sin j3) d cos (3 


J-[ 1-F(2*>] 


(84) 


where z = 2ka sin 1/2 9 . 


Lorentz-Lorenz Formula [8] 

The relation between polarizability and refractive index will now be derived. 

An external electric field stretches and orients molecular dipoles in a dielectric. 

The amount of polarization P is equal to the average electric dipole moment per unit 

volume. The polarization induces a field which modifies the relation between the intensity 
E and the displacement vector D. Hence, 
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(85) 


D - e 0 E + P 


The total polarization is caused by (1) polarization by stretching Pj and (2) polarization 
by dipole orientation P 2 . The polarization caused by elastically “stretched” molecules is 
given by 


P, = NaE' (86) 


where N = the number of molecules per unit volume, a = proportionality constant 
(polarizability), and E' is the effective internal electric field for the dipoles. 


The polarization caused by orienting the permanent dipole moments can be 
approximated by 


p 2 NE' 

P 2 = — — 

3 kT 


(87) 


where p is the molecular dipole moment, k is the Boltzmann constant, and T is the 
temperature. But for an oscillating field the relation becomes 



where r is the relaxation time for the dipoles, and r is related to the fractural torque on 
the molecules and the temperature. 

The Clausius-Mosetti theory gives the value of the internal field for a dipole as a 
function of the external field E and the polarization P. Consider a cavity within a 

dielectric. The local field in the cavity is made up of two parts: ( 1 ) the external field E 

and (2) the field produced by the polarization within the cavity. This second field can be 
calculated knowing that for a small area dA on the surface of a cavity the induced charge 
density is P cos 9 . For a spherical cavity the amount of charge induced in a spherical zone 

is 27rr 2 P sin d cos 6 d 9 , and the field at the center of the cavity is 1/2 e 0 P cos 2 9 sin 6 d 9 . 

/v 

Integrating over all the charges gives the total field as P/3 e 0 . Hence, the internal field is 
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E' = E + 


3 Co 


( 88 ) 


Using D = P + e 0 E = eE, the Clausius-Mosetti relation is obtained: 


E' = 


£ - Cq 


3 £o 



For polarization due to elastic stretching alone, 


(89) 


E 


— 3 go 


eo 


-5- + 2 

£o 


Na 


or in 


terms of the refractive index, m - y/e/eo. 



(Lorentz-Lorenz) 


Including permanent dipoles gives the following expression: 


(90) 


m 2 - 1 = _N_ a + J_ _El 1 

m 2 + 2 3 e 0 L 3 kf 1 + iur J 

showing the variation of the complex refractive index with temperature and frequence. 
For very high frequencies the second term in the bracket is negligible. Then the Rayleigh- 
Gans differential scattering cross section in terms of the refractive index is 


dn 

dO 



(1 + cos 2 6 ) R 2 (0,0) 


(91) 
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where 


Na = 3e ° \ 2 e 0 (m - 1) (92) 

m 2 + 2 m-> 1 

Scattering by Infinite Cylinders 

In conclusion we will consider the specific case of a thin infinite cylinder. According 
to Cooke and Kerker [9] , the radiation scattered from an infinite cylinder lies on the sur- 
face of a cone whose apical angle is twice the complement of the tilt angle. From Rayleigh- 
Gans scattering by cylinders (ka « 1 ), the expression for the scattered intensity is given by 
[7]: 


da 

dft 


= - 4 V g 2 - f 1 - )2 (1 + COS 2 9) ■ E 2 ^k£ sin cos (3 ^ 


(93) 


For infinite cylinders (£-><») E becomes significant only when cos (3 = 0, which includes 
9=0, the forward scattering direction. This yields by equation (77) 


cos 


a sin ^ 0^ = sin a cos ^ O^j cos <p 


(94) 


Using this restriction on a in the expression 


cos 7 = cos a cos 9 + sin a sin 9 cos <p 


or 


cosy = cos a ^2 cos 2 -y-l^ + sina ^2 sin —■ cos -y ^cos</> 


(95) 


where 7 is the angle between the cylinder axis and the scattering direction, gives the basic 
result that the scattered light is confined to the surface of a cone of apical angle 2a 
i.e., 7 = a. 
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Theoretically, one can show the formation of the scattering cone of the infinite 
cylinder by varying the length of a finite cylinder. Figure 9 shows the normalized scattering 
differential cross section 


do(0 ) 



dfi 



9 _ 

2 



(96) 



SCATTERING ANGLE 6 

Figure 9. The effect of varying the length £ of a finite cylinder 
of radius a = 0.03 /urn and refractive index m = 1.33 
( — £ = 0.1 /am, £ = 1 pm, £ = 10 /am). 

versus the scattering angle for the particular case where the scattered beams lie in the plane 
defined by the incident beam and the axis of the cylinder (Fig. 10). The cylinder of radius 
a = 0.3 /am has a tilt angle of 60 deg; hence, the scattered beams with scattering angles of 
60 deg and -120 deg coincide with the axis of the cylinder. For a length £ = 0.1 pm 
(£/a ~ 3), the peak due to the scattering cone has not appeared. A wide scattering cone 
peak is seen for the case £ = 1 pm (£/a ~ 33) at 0 deg and 120 deg. For £ = 10 pm 
(£/a ~ 333), the scattering pattern approaches that of an infinite cylinder with a scattering 
cone having an apical angle of 120 deg. The value of the normalized scattering differential 
cross section at d = 120 deg is nonzero, and not 1 as a result of a finite radius, i.e., 


1 

2 


F 2 


^ 2 ka sin 


^ (1 + cos 2 0) = 0.561 


The cone can easily be seen with a milliwatt helium-neon laser (0.6328 pm) and a 
slender wire [ 10] : Figure 1 1 shows the geometry of the experiment and photographs of the 
scattering pattern. Both the 25.4-pm-diam copper wire and the 1034-pm-diam wire shown 
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INCIDENT 

BEAM 


SCATTERED 



Figure 10. Cylinder scattering in the 
plane defined by the cylinder 
axis and incident beam. 


have a tilt angle of approximately 60 deg. 1 
The incident beam is made visible through the 
use of smoke. The scattering cone cross 
section is observed on a plane normal to the 
axis of the cylinder. Note that the incident 
beam lies on the cone in the forward scatter- 
ing direction. The thick wire is visible because 
of the light scattered by the smoke. The black 
patch seen in the photograph of the thin 
wire was used to reduce the amount of diffuse 
light coming from the screen. Cooke and 
Kerker [9] explain that the irregularities of 
the diffuse circles for the thick wire are 
caused by contamination on the wire, e.g., 
dust. The fine concentric circles which appear 
in the thin wire case are probably caused by 
irregularities on the surface of the wire. 


1 . Although Rayleigh-Gans theory is not applicable, the scattering geometry is equivalent 
to the Rayleigh-Gans case. 
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APPENDIX 


COMPUTER PROGRAM 


A listing of a computer program to calculate I, and I 2 for Ray leigh-Gans scattering 
for a cylinder and a disk is presented, followed by a sample output. The input/output 
information is given on the comment cards at the beginning of the program. A flow chart 
is presented in Figure A-l. For reference: 


R 2 (M) (A-l) 


I, = k 6 V 2 


( 


and 


I 2 = 

SD 

<*> 

O 

O 

w 

(A-2) 

I = 

(II +I 2 ) 
2 

(A-3) 

p = 

Ii -h 

(A-4) 


Ii - la 

da 

Ii + 1 2 

(A-5) 

dft 

2 k 2 
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SR UN » P CYLI NO. 4009 BP, . G AG ARYHST 55 1 . 3. 90 
3F0R.ISF MAIN 

£*****♦***«****************** 

C THIS PROGRAM COMMUTES II AND 12 FOR RAYlEIGH-GAN SCATTERING 

c sphere g=(3<stnu-ucosu>/u3> 

C OTSK Fr?6ESSlU/U 

C ROD E -S I NU /U 

C CYLINDER E*F 

C ELLIPSOI DS--USE G(U) WITH THE RADIUS REPLACED BY OC 

C (SEE VAN DE HULST PAGE 93 ) 

C I NPUT ***************** * 

C NTr N UMBER OF RUNS 

C INOEXlrlrSPHERE 

C =2:CYLINDER(DISK OR ROD AS LIMITS) 

C INDEX2=1- NO INTEGRATION OF BETA 

C :2-INTEGRATI0N OVER ORIENTATION 

C FM ^REFRACTIVE INDEX (REAL) 

C A rRAOIUS (MICRONS) 

C WAVE ^WAVELENGTH (MICRONS) 

C FLrLENGTH OF CYLINDER 

C BETArANGLE BETWEEN CYLINDER AXIS AND THE 

C BISECTRIX (OF SCAT DIRECT AND THE AXIS OF CYLINDER) 

C 0 U T°U T **************** * 

C ALR -ALPHA 

C XII =11 

C XI? r 12 

C P -POLARIZATT ON 

C THETA ^SCATTERING ANGLE (0-190 DEG. STEPR5) 

£**************************** 

14 F ORM AT ( 31 H * * * R A YLE I G H -G A N S SCATTERING***.//) 

5 F0RMATU8H REFRACTIVE IN OEX = . F7 . 3. / ) 

6 FORMAT ( 17H PARTICLE-RADIUS=.F7,3.8H lE NGT H= . F 7 . 3 . / ) 

7 FORMAT! 12H W A V EL ENGT Hr . F 7 , 3 . / ) 

8 F OR MAT ( 7H ALPHA-.F7.3./) 

9 FORMAT* 1H0.1 1H SCAT ANGLE. 7X.?0H INTENSITY FU NCT IONS . 1 9X , 

123H DEGREE OF POLARIZATION. //.11H T HE T A ( DE G ) . 7X . 3H I1.12X.3H 12. 
26X.8H 11*12/2. 13X. 2H P./) 

ID FOPMAT(lHl) 

11 F ORMAT ( 1H . I 5. 4X . 3F1 5.6. 4X , F9.6 ) 

2D FORMAT ( 8H INDEXlr .I3.22H 1=SPHERE 2 -CYlINDER) 

21 F ORMAT ( 8H T N OEX 2-. 13. 9H 1 z3ET A . F5 . 2 . 25 H 2- I NTEGR AT ION OVER BETA) 

1 FORMAT (13 ) 

22 FORMAT (213 » 2 F .1 0 . 6) 

READ (5 , 1 ) NI 

NNrl 

PI- 3. 1415 92 7 
PI2=2.*PI 

2 CONTINUE 

RE AD ( 5. 22) INDEX 1 . INDEX 2 . SET A . FL 
BETA-BETA *PI/18C. 

RE AO (5. 3 IF*. A. WAVE 

3 F0RMAT(3F10.3) 

FK-RI2/W AVE 

GO TO (23.24) .INDEX1 

23 V0L-(4.*PT*A**3)/3. 

GO TO 25 

24 V0L=(PI*A**2)*EL 

25 CONTINUE 
A3-1./3. 

C ALP FOR CYLINDER IS CALCULATED FOR EQUIVALENT SPHERE 
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ALP=PI2*( 3.*V0L/(4 • *P I > )**A3/WAVE 
CHI=(FK**6*V0L**2)*(( FM- 1 . > / ( 2 . *P I ) J * *2 
WRITE16.10) 

WRITE (6.14) 

WRITE(6.5 ) F M 
WRITE (6. 6) A. FL 
WRITE(6»7 ) WAVE 
WRITE (6* 8) ALP 
WRITE(6.20) INDEXI 
60 TO (26. 27 > .INDEX1 

27 WRITE(6»2I)INDEX2»BETA 
26 CONTINUE 

WRI T E ( 6 « 9 > 

THFTA -.0 
DO 4 1=1.37 

Z1 =2. *F«*A*ST N( THETA/2. ) 

Ul= Z 1*SIN (BETA ) 

Z2=FK*FL*5IN(THETA/2.) 

U2= Z 2*C OS (B ETA ) 

U=2.*ALP*STN( THETA/2.) 

GO TO (28 .29) .INDEXI 

28 ETA=(R(U) ) **2 
GO TO 30 

29 60 TO ( 31. 32 ) • I NDEX 2 

31 ETA = (F (U1 ) * E ( U2 ) ) * * 2 
60 TO 30 

32 CALL RBAR2(Z1.Z2.ETA> 

30 CONTINUE 

ZETA = (COS ( THETA ) ) * * 2 
X 1 1 = CHI * ET A 
XI2=XI1*ZETA 
XI 12=(XI 1+XI 2>/2. 

P=(XI1-XI2)/(XI1+XI2) 

NTH=5»I-5 

WRITE (6.1 1) NTH. XII . XI 2. XI 12. P 
THETA=THET A+5.*PI/180. 

4 CONTINUE 
NN = NN *1 

IF (NN-NI >2.2.13 
13 CONTINUE 
STOP 

FUNCTION G (U ) 

EPS= . 0D00D1 

IF (ABS(U)~E D S) 302.3C2.301 

302 G=l. 

GO TO 303 

301 G=( SIN ( U) -U*COS( U ) ) / (U**3 ) 

G = G*3 . 

303 CONTINUE 
RF TURN 

FUNCTION E(U') 

FPS=. 000001 

IF ( ABS(U)-EP5)312»312.311 
312 E=l. 

GO TO 313 
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311 E-SIN(U)/U 
313 RETURN 

FUNCTION F ( U ) 

FPS=. QOOOOl 

IF tABS<U)-E p S)322»322.321 

322 Frl. 

GO TO 323 
321 F-3SSL(U»1 ! 

F=F*2./U 

323 RFTURN 

FUNCTION BSSL(X.N) 

0-. 000001 
BJ-O. 

IF (N) 10. 20.20 
10 IER=1 
RE TURN 

20 IFfX >30.30.31 

30 IER=2 
RETURN 

31 IF (X- 15. >32. 32. 34 

32 NTEST-20.-.10.*X-X**2/3. 

GO TO 36 

34 NTESTr90.+X/2. 

36 IF (N-NTEST)40. 33.38 
38 IER- 4 
RETURN 
4.0 IERzO 
N1 -N + 1 
3PPE V- .0 

IF ( X- 5 . ) 50.60.60 
50 MArX+6. 

GO TO 70 

60 MA= 1 .4*X + 6D ./X 
70 MB=N* IFI X( X > /4+ 2 
MZERO=MAXP(«A.HB> 

MM AX^NTEST 

100 DO 130 M;M ZERO • MM A X • 3 
FM1-1 . 0 E - 2 8 
F Mz .0 
AL P H A - • G 

IF ( M-lM/2) *2 > 120 * 110 . 120 
110 jTr-1 

GO TO 130 
120 jTzl 
130 M 2- M -2 

00 160 K -1 . M2 
MK-M-K 

3MK=2.*F LO AT ( mk ) *FM1/X-FM 

FM=FM1 

FMlrBMK 

IF ( MK-N-l) 150 .140.150 
140 3J-SMK 
150 JT=-JT 
S=1+JT 

160 ALPHA- A L PH A ♦ 3 M K * S 
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BMK=2.*FM1/X-FM 
IF (N ) 180. 170. 1 30 
170 BJ=BKK 
180 AL p HA = ALPHA+'3MK 
BJ=B J/ALPHA 
BS SL = 8 J 

IF (ABS«B J-BPREV)-ABS(D*BJ) 1200.200. 190 
190 BPREV=BJ 
IEP=3 

200 RETURN 
END 

aFOR.IS R8AR2 

SUBROUTINE R3 AR2 C Z 1 . Z 2 . ET A 1 
C INTEGRATES OVER OR I EN AT I ON OF THE CYLINDERS 

DIMENSION Y ( IS 1 . A ( 1 6 1 
CALL INTIS(Y.A) 

ETA= .0 
DO 1 1=1.16 
F 1= ( Y ( I ) 1 

U1 =Z1 *S9 RT (1.-F1*F1) 

U 2= Z 2* F 1 

ET=(F { U 1 1 *E1U21 1 **2 
ETA=ET*A< I 1 .ETA 
1 CONTINUE 
RETURN 

FUNCTION G { U 1 
EPS= .000001 

IF (A8SlU)-E p S) 302.302. 301 

302 6=1. 

GO TO 303 

301 G= ( S IN ( U 1 -U *C0 S < U ) 1 / < U* * 3 1 
G = G*3 . 

303 CONTINUE 
RETURN 

FUNCTION E(U) 

E ° S= • 000 00 1 

IF (ABS(U)-E p S)312»312»311 

312 E = 1 • 

GO TO 313 
311 E = SIN ( U 1 /U 

313 RETURN 
FUNCTION F ( U 1 
E P S= .000001 

IF (ABS(U)-EPS) 322. 322. 321 

322 F= 1 . 

GO TO 323 
321 F=B S SL ( U » 1 1 
F =F *2 ./U 

323 CONTINUE 
RETURN 

FUNCTION BSSLtX.N) 

D=. 000001 
BU= 0 . 

IF (N) 10. 20.20 
10 IER=1 
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RETURN 

20 IF (X) iO. 30.31 

30 IER-2 
RETURN 

31 IF ( X -15. 1 32 2.34 

32 NTFS -20.+10. *X-X**2/3. 

GO T 36 

34 NTFS! =30 . ♦ X/ 2. 

36 IF(N-:!TEST) 4C.38.38 
38 IE R=4 
RETURN 
40 IE R=0 
N 1 = N + 1 
BPREV-.O 

IF( X-5. )5C. 60.60 
50 M A =X*6 • 

GO TO 70 

60 MA=1.4*X*60./X 
70 MB= N + IF I X IX 1/4+2 
MZERO=MAXQ<MA»MS> 

MMA XrNTEST 

100 00 130 M=MZERO. MMAX . 3 
FMl = 1 . 0E _ 28 
F M = .0 
ALPHAr.O 

IF (M-(M/ 21*2)120. 110. 120 
110 JTz-1 

GO TO 130 
120 JTr 1 
130 M2=M-2 

DO 160 Kr 1 « M2 
MK-M-K 

BMK=2. *FLOAT(MK ) *FM1/X-FM 

FM=FM1 

FMlrSMK 

IF (MK-N- 1 >150. 140.1 50 
140 BJ=BMK 
150 JT=-JT 
S = l + JT 

160 ALPHA=ALPH A+3MK *S 
8MK=2.*FM1/X-FM 
IF (N) 130. 170.130 
170 BJ=BMK 
130 AL°HA= ALPH A+ 3MK 
BJ-B J/A LPHA 
3SSL=3J 

IF ( A8S(B J-BPREV )-ABS«D*B J) 1200. 20 C. 130 
130 3°RE V =B J 
IER=3 

200 RETURN 
ENO 

3F0R.IS INT16 

SUBROUTINE INT1GCY.A) 

DIMENSION Y( 16 1 * A( 16 1 
C 16 POINT GAUSSIAti SUBROUTINE 
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Y AND A MUST BE DIMENSION IN CALLING PROGAM 

Y f 1 > — .95012510E-01 
Y ( 2 ) - 0.28160355 

Y ( 3) - 0.45301673 
Yf4 )- 0.61787624 
Y ( 5) - 0.75540441 
Y (6 ) - 0.86563120 
Y ( 7) - 0.94457502 

Y C 8 ) - 0.98940093 
A ( 1 ) - 0.18945061 
A ( 2 ) - 0.1 8260341 
MZ)- 0.1691 5652 
AC4)r 0.14959599 
A(5)r 0.12462897 
A(6)r . 35158512E-01 
A ( 7 ) - . 6 22 53 524E-0 1 
A(8)r .2715245 9E -01 
00 3 I - 9 f 1 6 

II- 17-1 

Y < I ) r - Y ( 1 1 > 

8 A (I )=A (II > 

RF TURN 
END 



♦♦♦RAYLEIGH-GANS SCATTERING*** 


REFRACTIVE INDEX: 1.300 


PARTICLE-RADIUS: .100 L ENGT H: 100 . 000 


WAVELENGTH: .500 


ALPHA: 11.417 

INDEX1: 2 1:SPHERE 2:CYLIN0ER 


T NnF X 2 : 1 

1 : 8 ET A • 1. 57 2: 

: TNT FGR AT TON OVER BETA 


SCAT ANGLE 

INTENSITY 

FUNCTIONS 


DEGREE OF P0IARIZATI0N 

THET A ( D EG ) 

LI 

12 

11 +12/2 

P 

0 

88601.623047 

88601.623047 

88601.623047 

.000000 

5 

88335.746094 

87664.736328 

38000.241211 

.003813 

LQ 

ft7*44. 11*7*4 


■ iniimiiB 

mm * r ~ ;l ' 

15 

86244. 449219 

80467 .165039 

83355.806641 

.034654 

20 

__ 8 44 65 r 4 75 586 


79525 - 182617 


25 

82245.773437 

67556.167969 

74900.970703 

.098060 


7 96 32 . 21 7773 

59724. 162598 


.142857 

35 

76678.210937 

51451 .852051 

64065.031 250 

.196881 



43097. 353516 

58269.522461 


45 

69983. 120117 

34991 .56QS47 

52487.340332 

.333333 

50 

66363.496094 

27419. 79&584 

46891 .646973 

.415252 

55 

62642.520996 

20608.758789 

41625.639648 

.504902 

GO 

5 88 76 * 954 5 90 

14719.239868 

36798 .097168 

.600000 

65 

55119.285645 

9844 .647095 

32481.966309 

.696920 


51416.681641 

6014.610291 

28715.645752 

i. 1 1 IU L f ^ ^ i ^ i; 

75 

47810.275391 

3202.682007 

25506.478516 

.874437 



■VCTPTHXIFa 


.941458 

85 

41018.243164 

311.579933 

20664.911377 

.984922 


■ mmiiiiiw 



1 .DOODOO 

95 

34943.179687 

265.432404 

17604.305908 

.984922 

R£S5ii^HT«T«^BSEEi 


^Eril^rwiv 

■ IB 1 1 U ■ 


105 

29689.675781 

1988.830307 

15839.252930 

.874437 

110 

27381.871826 


15292.470581 

.790546 

115 

25284.722656 

4516.006653 

14900.364624 

.696920 

Ml liM 111 'I' " 1" 

■ZuTICTXTiTT^B 

■lTTTTWTTTTTI 



125 

21700.095459 

7139.110352 

14419.602905 

.5049D2 


■ Mimimiw 

RlfllKViUM 


.415252 

135 

18874.645996 

9437.319946 

14155.982910 

.333333 

140 

17723.308105 

10400.461060 

14061.884521 

.260379 

145 

16733.472900 

1 1228 .325806 

13980.899292 

.196881 

150 

1 58 86 .047119 

11922.032349 

13909.039673 

.142857 

155 

15202.678955 

12487.383545 

1 3845.031250 

.098060 

160 

1 4645.975708 

12932.719727 

13789 .347656 

.062122 

165 

14219.686890 

1 3267 .146606 

13743.416748 

.034654 

170 

13918.827026 

13499.121460 

13708 .974243 

.015308 

175 

13739.775757 

13635.405762 

1 3687.590698 

.003813 

U51L 

13680.336426 

13680.336426 

I368D .336426 

.000000 
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FUNCTIONS 

G(u) = 3 (SINu~uCOSu)/u 
F(u) = 2J i (u)/u 
E(u) = SINu/u 

Jj (u) = BESSEL FUNCTIONS 



Figure A-l. Flow chart for Rayleigh-Gans scattering program. 
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